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ABSTRACT 



We used FUSE to observe ultraviolet emission from diffuse VI in the hot gas 
in the Galactic halo. By comparing our result with another, nearby observation 
blocked by an opaque cloud at a distance of 230 pc, we could subtract off the 
contribution from the Local Bubble, leading to an apparent halo intensity of Iovi 
= 4680l^g photons cm 2 s 1 sr 1 . A correction for foreground extinction leads to 
an intrinsic intensity that could be as much as twice this value. Assuming T ~ 3 x 
10 5 K, we conclude that the electron density, n e , is 0.01 — 0.02 cm -3 , the thermal 
pressure, pj k, is 7000—10, 000 cm -3 K, and that the hot gas is spread over a length 
of 50-70 pc, implying a small filling factor for O VI- rich gas. ROSAT observations 
of emission at 1/4 keV in the same direction indicate that the X-rays are weaker 
by a factor of 1.1 to 4.7, depending on the foreground extinction. Simulated 
supernova remnants evolving in low density gas have similar O VI to X-ray ratios 
when the remnant plasma is approaching collisional ioinizational equilibrium and 
the physical structures are approaching dynamical "middle age". Alternatively, 
the plasma can be described by a temperature power-law. Assuming that the 
material is approximately isobaric and the length scales according to T@d In T, 
we find (3 — 1.5 ± 0.6 and an upper temperature cutoff of io 6 - 6 (+°- 3 >-°- 2 ) K. The 
radiative cooling rate for the hot gas, including that which is too hot to hold O VI, 
is 6 x 10 38 erg s _1 kpc~ 2 . This rate implies that ~70% of the energy produced in 
the disk and halo by SN and pre-SN winds is radiated by the hot gas in the halo. 



-2- 



Subject headings: Galaxy: general — Galaxy: halo — ISM: general — ultraviolet: 
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Introduction 



Ultraviolet and soft X-ray observations by Copernicus, Voyager, HUT, ORFEUS, FUSE, 
SPEAR, the Wisconsin All Sky Survey, EE AO, ROSAT, Chandra, and XMM detected O VI 
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and 3/4 keV soft X-rays trace 10 6 to 10 7 K gas (IMazzotta et al.lll998l ; iRaymond fc Smith 



19771 ). It is natural to anticipate that the O Vl-rich and 1/4 keV emitting plasmas are 
causally, and perhaps physically, associated because million degree gas eventually cools into 
the 300,000 K regime. As it cools, its O VII ions recombine with electrons to make O VI ions. 
Furthermore, as time goes by, interfaces between 10 6 K gas and warm ionized gas (~ 10 K) 
should develop at temperatures that are between these two values and exhibit intermediate 
levels of ionization. Thus, it is reasonable to expect the 1/4 keV emitting gas to be clad 
in an O Vl-rich sheath, irrespective of whether the hot extraplanar gas is due to fountains, 
superbubbles, or extraplanar super nova (SN) explosions. Su c h bub ble-sheath geometries are 
present in the disk simulations of lAvillez fc Breitschwerdtl (120041 ). as are T ~ 3 x 10 5 K 
bubbles without T ~ 10 6 K cores. 

In this paper, we focus on hot gas above the Galactic disk. Even though the nearest 
part of this region may be only a few hundred parsecs above the midplane, we will follow the 
X-ray community's tradition by referring to this region as the halo. There are a number of 
fundamental issues that we plan to address about the character of the material that creates 
O VI and X-ray emissions: 1.) Are the source structures in the halo young, middle aged, or 
ancient? 2.) What is the temperature distribution function of this gas? and 3.) How long 
is this plasma's cooling timescale and how does the energy loss rate compare with the rate 
of energy injection from supernova explosions in both the disk and halo? 

Our analyses are based on sets of Far Ultraviolet Spectroscopic Explorer (FUSE) O VI 
and Roentgen Satellite (ROSAT) 1/4 keV emission data for I = 278.7°, b = —47.1° and 
/ = 278.6°, b = —45.3°. We focus much of our attention on the O VI emission data, which 
can be compared to the soft X-ray emission to synthesize a long baseline spectrum. In turn, 
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this spectrum defines the shape of the plasma's temperature distribution function and can 
be compared with models of supernova remnant (SNR) evolution to constrain the age of a 
remnant, on the assumption that the hot gas in our line of sight is dominated by the products 
of one such remnant. 

We have no direct measurement of the O VI column density along our line of sight. 
Nevertheless, observations of O VI absorption features toward bright extragalactic sources 
at only moderate angular separations from our direction give an approximate guide on the 
value of the O VI column density, iVovi- This information is useful, since both O VI and soft 
X-ray emission intensities are equal to (1/4-71") J n e rii < av > (T)dl, where n e is the local 
electron density, rii is the local density of ions, and < av > (T) is a temperature-dependent 
emissivity function. In contrast, the O VI column density is simply equal to the line of 
sight integral of the ion density. The column density integral is not additionally weighted by 
electron density, and its temperature dependence is simply that due to the changes in the 
ion fraction. As a result, the ratio of the O VI intensity to column density can be used to 
determine the electron density and from that, the thermal pressure. 

Every observation of the halo's intensity is contaminated by photons produce d elsewhere 



along the line of sight. The most obvious non-halo sources are the heliosphere (ILallement 



20041 ) and the Local Bubble (also called the Local Hot Bubble), where the latter is an irreg- 



ularly shaped pocket of X-ray emissive, interstellar g as surrounding the Solar neighborhood 



out to a distance of ~100 pc (I Cox &: Reynolds! 1 19871 ). The combined Local Bubble and he- 



liospheric emission must be subtracted from that found for observations of the high latitude 
sky in order to determine the halo's intensity. To accomplish this, we selected two adjacent 
lines of sight, one that is nearly free of foreground obscuration and another that has a cloud 
of neutral gas, which blocks the more distant halo emission but not the foreground emission. 
By evaluating the difference between the two intensities, we can measure the emission that 
arises only from the halo. This procedure is called the "shadowing technique." 

For our shadowed line of sight, we used the strength of O VI emission found along 
the path to a nearby , opaque filament. The filament is located 230 ± 30 pc from the Sun 



( IPenprase et al.lll998l ). so is positioned beyond, but not far beyond the Local Bubble. FUSE 
has already observed a dense part of the filament at I = 278.6°, b = —45.3°. The observation 
yielded a 2a upper limit on the O VI doublet's in tensity, including the systematic uncertain- 



ties, of 800 photons s 1 cm 2 sr 1 (jSheltonl|2003l ). Ignoring the systematic uncertainties and 



using the la statistical uncertainties as error bars leads to a best fit value of 30 ±340 photons 
s _1 cm -2 sr -1 (as determined by one of us (RLS) while preparing the original article). To 
measure the combined halo and Local Bubble O VI intensity for this paper, we used FUSE 
to observe an unobscured line of sight only 2° from the filament {I = 278.7°, b = —47.1°; 
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see Figured]). In order to estimate the halo's soft X-ray count rate, we perform a similar 
shadowing analysis using ROSAT data for the same pointing directions. We also subtract 
the extragalactic soft X-ray contribution and account for the variation in optical depth with 
photon energy. By carefully isolating the halo emission, we expect to circumvent the amal- 
gamation problem seen in plots of total O VI intensity versus total soft X-ray countrate. 

In Section [2] of this paper we present our measurements of the VI emission produced 
along the unobscured, off- filament line of sight. In Subsection 12. 4[ we subtract the local 
O VI intensity from our measurement of the unobscured line of sight, yielding the intensity 
originating in the halo. We determine the halo's O VI column density from published data 
for nearby directions in Subsection 12.51 In Subsection 12.61 we estimate of the halo's 1/4 keV 
surface brightness by performing a shadowing analysis with ROSAT data from the filament 
and off-filament directions. We then use the halo's O VI intensity and column density 
determinations and treat the O Vl-rich gas as if it is isothermal in order to estimate the 
density and pressure of the emitting material in Subsection 13.21 We compare the O VI and 
1/4 keV results with simulation predictions in order to estimate the plasma's maturity in 
Subsections 13.41 We compare the O VI and 1/4 keV results with analytic functions in order 
to estimate the plasma's temperature distribution function, pressure, and cooling rate in 
Subsection 13.51 The results are summarized in Section HI 

2. Observations, Data Reduction, Spectral Analysis, and Determination of 

the Halo's O VI Intensity 

2.1. Observations 

FUSE observed the blank sky toward I = 278.69°, b = -47.08° (R.A. = 03 h 20 m 10 s , decl. 
= — 62°26'29") once in July and twice in November of 2002, for Guest Observer program 
C153. These observations exposed the LiF 1A detector segment, the segment used for the 
subsequent O VI signal extraction, for 53, 104, and 61 ksec, respectively. Of these durations, 
6, 61, and 26 ksec, respectively, occurred while the satellite was in the night portion of 
its orbit. FUSE also made a non-proprietary "Early Release Observation" of this direction 
in December 1999 and archived the data under program identification number X0270301. 
However, detector 1 was de-powered during most of the observation time. 
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2.2. Data Reduction 

During most of the Guest Observer exposures, data were taken on all four detector 
segments (1A, IB, 2 A, and 2B). The data from each detector segment were taken in time- 
tag mode where the photoevents from each detector are listed in the sequence that they 
were recorded, with their times and locations on the detectors where they were found. We 
concatenated the photon lists obtained from the individual exposures, creating lists for each 
detector segment and observation. 

We then ran the concatenated photon lists thro ugh version 3.1.1 of the CalFUSE pipeline 



( jSahnow et al.ll2000l : iDixon. Kruk fc Murphyil2002l ). with corrections enabling us to use the 



extended source extraction window. We did not use automatic background subtraction. 
For the LiF 1A extractions, we set the lower and upper pulse height limits to 4 and 25, 
respectively. For the other detector segments, we used the default pulse height limits. This 
process yielded several spectra, because in addition to recording emission with four detector 
segments, FUSE records light reflected off of each of two types of coatings (lithium fluoride 
(henceforth called LiF) and silicon carbide (henceforth called SiC)) and four gratings. The 
multiplicity enables the telescope to see eight overlapping wavelength regimes with various 
efficiencies. Furthermore, each mirror's focal plane assembly has four apertures. Thus, each 
detector segment is able to receive light from each of the four apertures simultaneously and 
without overlap. The spectra observed through the largest aperture (LWRS: 30" x 30") in 
the LiFlA channel provides the greatest signal to noise and the greatest effective area in the 
1030 to 1040 A regime. We will use these spectra for the following O VI data reduction, and 
will use the spectra from the other channels when searching for other cosmic emission lines 
outside of the 1030 to 1040 A bandpass. Our avoid ance of the SiC channels overcome s the 



solar O VI line contamination problem discussed by iLecavelier des Etangs et al.l (120041 ) 



The wavelength scales for each spectrum require re-calibration, which we determined by 
comparing the reported heliocentric wavelengths of the observed geo c oronal emission lines 



around 1027 and 1040 A with their rest wavelengths (jMortonl Il99ll ) converted from the 
geocentric to the heliocentric reference frame. We fit a linear function of wavelength to the 
difference between the expected and observed heliocentric wavelengths, then subtracted this 
function from the observed wavelengths. We then combined like spectra and shifted to the 
Local Standard of Rest (LSR) reference frame. 
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2.3. Spectral Analysis 

Figure [2] displays the satellite-night and day+night LiF 1A spectra for wavelengths 
between 1027 and 1045 A. The spectra reveal O VI emission in the 1032 and 1038 A lines 
(rest wavelengths: 1031.93 and 1037.62 A) and C II* emission at 1037 A (rest wavelength: 
1037.02 A). An additional emission feature appears in the day+night spectrum at 1031 A, 
but does not appear in the night-only spectrum. It also appears in other long-duration 
FUSE daytime spectra but not in their nighttime counterparts. We suspect that the 1031 A 
feature is the second order diffraction peak associated with atmospheric He I, whose rest 
wavelength is 515.62 A. Other He I second order lines in our bandpass include 2 x 522.21 A 
and 2 x 537.03 A, both of which appear in our daytime spectrum but not our nighttime 
spectrum. Along with the emission lines is also a continuum, which arises from radioactive 
decay within the detector, high energy particles, and scattered terrestrial airglow photons. 
The last of these is strongest during the daytime portion of the orbit. 

We used two methods to measure the stre ngths of the i n terste llar emission features. 



With the first method, which parallels that of IShelton et all (120011 ). we fit a second or- 
der polynomial to the continuum of the counts spectrum. (Note that "counts" are called 
"weights" in CALFUSE version 3.1 output files.) We then subtracted the continuum from 
the spectrum and searched for wavelength regions with large numbers of counts. Due to 
continuum subtraction and statistical variation, a few of the pixels within weak emission 
features have negative numbers of counts and some pixels outside of the emission features 
have positive numbers of counts. Although this noise complicated our effort to determine 
the boundaries of emission features, we set the boundaries of the emission features so as 
to exclude neighboring regions exhibiting near random fluctuations. We then summed the 
residual counts in each emission feature, divided by the instrumental effective area, and in- 
corporated geometric factors in order to convert the measurements to units of intensity. Each 
measurement's statistical uncertainty was calculated from the square root of the number of 
counts in the original spectra between the upper and lower wavelength range of the feature. 

In the day+night spectrum, the O VI 1032 A emission line and the daytime 1031 A 
feature are too close together to be separated with this algorithm. Therefore, when we 
analyzed the 1032 A emission line in the day+night data, we used the wavelength range found 
from the nighttime data (1031.64 A to 1032.45 A). As can be seen in Figure [U the lower end 
of this range overlaps very little with the He I feature. Furthermore, the measured intensity 
in the day+night 1032 A O VI feature is less than that of the night time feature, leading us to 
believe that little or no He I emission was attributed to the day+night 1032 A feature. Tabled] 
lists the measured intensities and 1 sigma statistical uncertainties of the O VI 1032 and 1038 
A lines, the C II* line, and the anomalous 1031 A feature. For presentation purposes, we 
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round the intensities as well as the Local Bubble intensity in §2.41 to 10s of photons sec _1 
cm -2 sr -1 . We use the unrounded values when calculating resulting quantities. Given the 
uncertainties in ap erture size and flux calibrations, additional systematic uncertainties of 
14% are expected (IShelton et all 120011 ). Placement of the continuum is another source of 
uncertainty. Since the continuum was independently determined during each analysis, the 
variation within our set of measurements provides an estimate of the size of this uncertainty. 
The standard deviation among the 4 measurements of the 1032 A emission line (Table [T|) is 
220 photons s^ 1 cm -2 sr -1 while the standard deviation for the 1038 A line is 470 photons 



s 1 cm 2 sr 



In order to determine the central wavelengths of our irregular, weak emission features, 
we fit each residual spectrum with a Gaussian and a quadratic. Generally, the fitting routine 
modeled the emission line with the Gaussian and the residual background with the quadratic. 
Occasionally, as occurred in the O VI 1032 A fits, the Gaussian extended beyond the bulk of 
the emission line in order to include nearby weak residual background emission. Thus, the 
reported wavelengths for the day+night O VI 1032 A Gaussian fit and the nighttime O VI 
1038 A Gaussian fit are higher and lower, respectively, than one would estimate by eye. See 
Table [2] for the tabulated results and note that the corrected FUSE scale is accurate to only 
~ 10 km sec -1 . 



In the second method, which parallels that of iDixon et al.l (120011 ). we fit the intensity 
spectrum's continuum with a straight line and its emission features with Gaussian shapes 
that had been convolved with a 106 km sec -1 top-hat function (simulating the finite width of 
the LWRS aperture). Because the O VI 1032 A emission line and the daytime 1031 A emission 
feature are so close together in the day+night spectrum, they were fit simultaneously. The 
O VI 1038 A emission line and the C II* emission line required simultaneous fitting for the 
same reason. The fits occasionally included outlying emission, which shifted the reported 
centroid wavelengths and increased the areas unde r the curves. The la error bars for each 
model parameter were determined as described in IDixon et al.l (120011 ). See Tables [1] and [2] 
for the measured intensities and velocities. 

The O VI and C II* appear to be moving very slowly relative to the Local Standard 
of Rest. However, the 1032 A signal is offset by about 30 km s _1 from the 1038 A signal, 
even in the fits that were not corrupted by outlying emission, the night-only 1032 A and 
the day+night 1038 A fits. The difference may be due to inaccuracies in the wavelength 
scale, being as the observed difference is somewhat more than twice the uncertainty in the 
velocity of each feature. It could also be due to variations in self- absorption optical depth 
as a function of wavelength. Larg er differences have b een seen in other spectra, for example 
the Virgo and Coma spectra of IDixon et al. (hooih . In the night-only 1032 A and the 
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day+night 1038 A fits, the o of the cosmic emission lines' Gaussian profiles were 47 and 
39 km s _1 , respectively. The resulting full widths at half maximum are 111 and 92 km s _1 , 
respectively, while the Doppler parameters are 65 and 55 km s _1 , respectively The true 
widths of the cosmic emission lines should be slightly narrower than these values because 
the fitting method accounts for only the broadening due to a finite slit width and not any 
other form of instrumental broadening. 

As a result of applying two analysis methods to two dataset variants (the day + night 
data and the night only data), we have 4 measurements for each of the O VI 1032 A, O VI 
1038 A, and C II* 1037 A emission lines. For the most part, the la measurements overlap. 
However, for each cosmic line, 3 of the 4 intensity measurements are similar to each other 
while the 4th measurement is an outlier. From the 3 clustered and unrounded measurements, 
we calculate average intensities of the O VI 1032 A line, 1038 A line, and the C II* line. They 
are 3190±450, 1520±350, and 1600±330 photons sec -1 cm -2 sr _1 , respectively. 

Our O VI signals are neither especially bright nor especially dim when c omp ared with the 



mid la titude (35° < b < 60°) FUSE detections in the lOtte fc Dixonl (120061) a nd lDixon et al. 



(120061 ) survey. The larger of these compilations, that of iDixon et al.l (120061 ). lists 17 sight 
lines having 2 o and 3 a detections of moderate velocity (-100 km s" 1 < v < 100 km s _1 ) 
O VI 1032 A emission. Of these 17 signals, 10 were brighter and 7 were dimmer than our 
1032 A signal. 

We also examined the data for signs of other cosmic emission lines. No significant signals 
were found. Our search for C III emission at 977.02 A in the satellite-night SiC IB spectrum 
yielded an intensity ±1(7 of 1300 ± 1740 photons sec -1 cm -2 sr _1 . Similarly, our search for 
N II emission at 916.71 A in the same spectrum yielded 510 ± 1250 photons sec -1 cm -2 
sr _1 . Although neither species is observed in our dataset, both are observed at the > 2cr 
level in the on-filament observation (C III: 4700 ± 1300 and 2600 ± 1000 photons sec -1 cm -2 
sr _1 in the nighttime SiC IB and SiC 2A spectra, r espective l y; N II: 2300 ± 1100 photons 



sec -1 cm~ 2 sr _1 in the night time SiC IB spectrum. ISheltonl (120031 )). Perhaps the observed 
photons came from the filament itself. Though, if this is the case, then the C III in tensity 
is rem arkably bright, considering that theoretical predictions for evaporating clouds (jsiavin 



19891 ). predict an order of magnitude dimmer emission. 



2.4. Determination of the O VI Intensity from the Halo 

Along our line of sight, the combined intensity arising from both members of the O VI 
doublet is 4710 ± 570 photons cm~ 2 s _1 sr _1 . The observed intensity can be apportioned 
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between the Local Bubble (plus heliospheric) and 'halo' components if we assume that the 
local emission is approximately constant over small angular separations. The observation 
of an opaque filament o n a nearby di rection, £ = 278.6°, b = —45.3°, set upper limits on 
the local O VI intensity (jSheltonl 120031 ) . For the local contribution, we take the tightest 1 a 
upper limit (due to random variation only) on the doublet from the day+night data for the 
1032 A emission line and the assumption that the 1038 A line is half as strong as the 1032 
A line. This doublet intensity is 30 ± 340 photons cm -2 s _1 sr _1 . However, since the Local 
Bubble cannot produce negative photons, we adjust the estimated intensity arising within 
230 ± 30 pc of the Earth (thus within 160 ± 20 pc of the Galactic plane) to 30^3q° photons 
cm -2 s" 1 sr _1 . Thus, for our sight line, the remainder and vast majority of the observed 
intensity (4680^ggQ photons cm -2 s -1 sr _1 ) is produced in the halo. 



Note that the filament data were later included in the iDixon et al.l (120061 ) catalog. To 
understand the lDixon et al.l (120061 ) reanalysis, the reader needs to first know that the filament 
was observed in 5 sets of grouped LWRS exposures, taken over a 2 year period, towards 3 very 
similar directions on the sky {{£ = 278 58°, b = -45.3 1°), {£ = 278.59°, b = -45.30°), and 
(£ = 278.63°, b = -45.31°)). Figure 3 of lsheltonl d2003h shows that a ll three directions point 
toward a dense knot in an infrared-bright filament. ISheftonl (120031 ) separately plotted the 
spectra for each of the three directions, found no important systematic differences between 
the spectra, then added the spectra and searched for emission in the O VI 1032 and 1038 A 
lines. No emission was found near the Milky Way's Local Standard of Rest (LSR) velocity 
in either line, in either the data taken during the satellite night portion of the orbit or 
data taken during the the satellite day + night portions of the orbit. In our shadowing 
analysis (previous paragraph), we have adopt ed the tightest la upper limit for the doublet, 
excluding systematic uncertainties, from the ISheltonl (120031 ) analysis. It was derived from 
the day+night data for the O VI 1032 A region. 



The later analysis by IDixon et al.l (120061 ) differed from that of ISheltonl (120031 ) because 
it included only data taken during orbital night, searched for only the 1032 A emission 
line, and used a different version of the pipeline, different choices of pulse height cutoffs, 
and a different spectral fitting algorithm. They noticed a high velocity feature in their 
combined dataset, but determined that it was primarily associated with their data taken 
towards £ = 278.59°, b = -45.30°. In that direction, their O VI 1032 A feature's velocity 
and intensity are 206 ± 13 km s -1 , and 3.0 ± 0.6 x 10 3 photons cm -2 s~ x sr _1 , respectively. 
In contrast, neither of their spect ra in the other dir ections revealed statistically significant 
O VI 1032 A emission. From this. IDixon et al.l ( 20061 ) concluded that O VI in the Magellanic 
Stream moving at ~ 200 km s _1 relative to the LSR had been viewed through a previously 
unknown hole i n the filamen t and that the other observations viewed opaque portions of the 
filament. When lDixon et al.l (120061 ) excluded the red-shifted feature from their analysis of the 
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combined dataset, they derived a 2a upper limit on the O VI 1032 A emission of 600 photons 
cm -2 s _1 sr _1 (see the last par agraph of thei r Appendix A. 2 and note that their naming 
convention differs from that of Sheltonl (j2003 )). which is roughly similar to the 2er upper 
limit reported by ISheftonl (120031 ). Giv en that they d id not find LSR rest frame emission, 
they did not modify the conclusions of ISheltonl (120031 ) regarding the physical conditions in 
the Local Bubble. 



If the high velocity O VI reported by lDixon et al.l (120061 ) were to lie along the off-filament 
line of sight analyzed in this paper, then its 1032 A emission would be 90% as bright as and 
0.6 A longwards of the O VI 1032 A feature we observed. As can be seen in Fig. [2J our 
spectra are dim around A = 1032.64 A, the wavelength of O VI traveling at 206 km s _1 . Our 
wavelength calibration in this part of the spectrum is good to within ~0.034 A (~10 km s -1 ); 
therefore, the observed velocity o f the O VI 1032 A f eature in our spectrum is inconsistent 
with the reported velocity of the iDixon et al.l (120061 ) fe ature. Thus , we c onclude that the 
high velocity, presumably extra-galactic feature seen by IDixon et al.l (120061 ) in the 12050501 
+ 12050510 data has not appeared in the off-filament data and does not affect our analysis 
of the halo's emission. 

Next we consider the effects of extinction. Our off-filament line of sight is only mildly 
extincted, but we obtain so mewhat different results from differen t datasets. According to 
the DIRBE and IRAS data (jSchlegel. Finkbeiner. &: David Il998l ). its color excess is E(B- 



V) = 0.0217, implying that N(R I) = 1.06 x 1 °cm if we u s e the empirical relation 
A(H l)/E{B — V) — 4.93 x 10 21 cm" 2 derived bv biplas fc Savage! Jl994h . Thus, we expect 



that 25% of the photons originating beyond the obscuring material have been extincted 
(|Fitzpatrick!ll999l) . However, according to the Leiden- Argentine-Bonn Survey (IKalberla et al 
20051 ). the column density of neutral hydrogen, Ahi, is 2.04 x 10 20 cm -2 on the two nearest 



l ines of sight ((I, b) = 27 8.69, -47.08 and 278.80, -47.0 0), equating to a c olor excess of 0.0414 



(IDiplas fc Savage 



19941 ) and an extinction loss of 42% (jFitzpatrick!ll999l ) . We recognize that 



this may slightly overestimate the amount of hydrogen along our line of sight because some 
small portion of the 0?6 (HPBW) diameter beam could be responding to the emission from 
the filament. We will take Ahi = 2.0 x 10 20 cm -2 as the upper end of the possible range. 

For a lower limit to the absorption, we consider the extreme case that a large fraction 
of the neutral hydrogen might be beyond the O Vl-emitting gas (almost all of which must 
be more distant than the shadowing filament). A limit to that fr action is defined by an 
estimate for the amount of H I that exists within the Local Bubble. lLallement et al.l (120031 ) 
found that well in front of the distance to the filament (230 ± 30 pc) the equivalent widths 
of intervening Na I D absorption in the spectra of stars in the general vicinity of our sight 
line are at least 20 mA and perhaps even as large as 50 mA (corresponding to Ahi = 0.5 to 
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2.0 x 10 20 cm -2 , respectively.) Thus, there must be some extinction between our location and 
the O VI ions. Allowing for small scale variations in ISM column densities and for possible 
revisions in the distance estimates, we take Ahi = 0-5 x 10 20 cm -2 as our lower limit. Also, 
we take Nm = 1.0 x 10 20 cm -2 as our nominal value. For Nm = 1.0 x 10 20 cm -2 , the intrinsic 
intensity, Jovi, is 6110+ggo (11. x 10~ 8 erg cm -2 s -1 sr _1 ). For the extreme range of 
possible extincting column densities, Nm = 0.5 to 2.0 x 10 20 cm -2 , the intrinsic intensity, 
lovi, is 53501^8 photons cm" 2 s" 1 sr" 1 (10.3±i;§ x lO" 8 erg cm" 2 s" 1 sr" 1 ) to 7960^ 
photons cm -2 s _1 sr _1 (15.3l 2 *2 x 10 -8 erg cm" 2 s _1 sr _1 ). 



2.5. An Estimate for the Halo's O VI Column Density 



In order to determine a characteristic electron density for the hot gas, we must know 
the halo's O VI column density (iVovi)- To estimate Aovi along our line of s ight, we have 
at ou r disposal only determinations made in other directions about 15° away (IWakker et al 
20031 ). We make use of 4 such measurements and evaluate their average column density. 
This average equals 2.34 x 10 14 cm -2 towards (1) NGC 1705, (2) Fairall 9, (3) several targets 
within the LMC, which we treat as a single data p oint, and (4) severa l targets within the 
SMC, which we also treat as a single data point (jSavage et al.ll2003l ). While this mean 
value represents our best estimate for our direction, we recognize that O VI column densities 
vary markedly over angular separations of 15°, and thus our adopted value could deviate 
substantially from the true value. 

We now estima te the size of the error in our determination. By examining Figure 11 of 
Savage et al.l (120031 ) . we can better estimate the expected deviation from the parent popula- 
tion than simply evaluating the internal dispersion of the nearest 4 lines of sight considered 
here. The average difference i n column den s ity of individual sight lines in pairs with ~ 15° 
separations in the data set of ISavage et al.l (120031 ) is 35% of the mean of each pair. Thus, 
on average, each determination differs from the mean of the two by 17.5%. The root-mean- 
square (rms) deviations should be \/tt/2 times this value if the distribution is Gaussian, 
and individual samples drawn from the parent population have dispersions that are a factor 
of y/2 larger, leading to a fundamental rms deviation of 31% from what we could consider 
to be a "true" overall mean value over some arbitrarily large sector of the sky. The mean 
of 4 randomly drawn samples from the parent population should deviate from the overall 
mean by only half as much (15.5%). Thus, our estimate for iVovi along our line of sight 
based on the other 4 measurements could be in error by an amount equal to the error in 
the mean (15.5%) combined in quadrature with the intrinsic fluctuations in the individual 
lines of sight (31%), giving a total uncertainty of 34.7%. Thus, we estimate the O VI column 
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density and 1 a error for our line of sight should equal 2.34 ± 0.81 x 10 cm 



From this value, we must subtract the Local Bubble's column density, which we estimate 
from i ts O VI volume density and radius. Three of the white dwarfs in the ISavage fe Lehner 
(120061 ) survey of nearby stars are within 30° of our pointing direction. The volume densities 
on these sh ort and nearby lines of s ight range from 2.01 to 10.6 x 10 -8 O VI ions cm -3 . From 
Figure 6 in lLallement et al.l (120031 ). we take the radius of the Local Bubble in our direction 
to be 80 to 140 pc. From the minimum radius and O VI volume density, we estimate the 
minimum Local Bubble O VI column density to be 5.0 x 10 12 cm -2 , and from the maximum 
radius and O VI volume density, we estimate it to be 4.58 x 10 13 cm -2 . Conservatively 
adopting these values as the Local Bubble's average column density ± 1 a and subtracting 
them from the column density toward extragalactic targets yields the halo column density, 
iVov! = 2.09 ± 0.84 x 10 14 cm" 2 



2.6. An Estimate for the Halo's X-Ray Brightness 

We now use results from the ROSAT All Sky Survey to find X-ray intensities that ap- 
ply to our two lines of sight. We extract the 1/4 keV X-ray count rate from the survey 
data for the Rl and R2 bands (~ 110 to 284 eV and ~ 140 to 284 eV, respectively). 
A 0.4° radius disk centered on our direction I = 278.7°, b = —47.1°, has a combined 



ROSAT Rl + R2 countrate of 1322 ± 55 xlO -6 (ISnowden et all (119971 ). retrieved with 



the X-Ray Background Tool: |http : / / heasar c . gsf c . nasa . gov / cgi-bin / Tools /xra ybf / xr aybg. pi ) . 



This countrate includes contributions from the halo, Local Bubble, heliospheric charge ex- 
change, and extragalactic objects. Assuming that the Local Bubble contribution is roughly 
uniform over small angular scales and that the heliospheric contribution is roughly time 
invariant, we take the Local Bubble plus heliospheric contributions to equal the coun- 
trate seen toward the nearby filament (534 ± 131 xlO -6 counts s -1 arcmin -2 for a 0.1° 
disk centered on I = 278.6°, b = —45.3°). In order to estimate the extragalactic contri- 
bution, we draw upon the parameterizations of Miyaji et al. (1998). They found that 
1(E) = 10.0 ±0.5 x E- 1A2±om and 1(E) = 11.5 ±0.7 x £-^±0.03 photons cm -2 k -i ke y-i 
gr -i f or pi 10 t on energies, E, between 0.1 and 10 keV for the two fields they studied. We 
take the average: 1(E) = 10.75 x E~ 1A2 photons cm -2 k -1 keV -1 sr -1 . When absorbed by 
galactic material (here we take the total column to be as high as Nm = 2.0 x 10 20 cm -2 ) and 
convolved with the ROSAT response function, such a spectrum yields 55.6 xlO -6 counts s -1 
arcmin -2 in the R1±R2 bandpass, hereafter termed the R12 bandpass. Subtracting the local 
and extragalactic contributions from the observed countrate along our line of sight yields 
a halo countrate of 732 ± 142 xlO -6 counts s -1 arcmin -2 . This rate applies if emission 
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originates below all of the observed neutral hydrogen. However, we believe that the halo 
emission originates above an absorbing column of 0.5 to 2.0 x 10 20 cm -2 . In this case, the 
intrinsic countrate would be 779 ± 151 to 4770 ± 930 counts s _1 arcmin -2 . 



For the analysis that will be presented in § I3.5[ we will also require the halo's R6 + 
R7 band (hereafter termed the R67 band) countrate . The filament is too opti c ally t hin at 



1.5 keV for a useful shadowing analysis. However, iHenley. Shelton fc Kuntzl (120061 ) have 
performed spectral fits with the RO SAT and XMM data for the on-filament and off-filament 
directions. Their analysis takes into account the extragalactic and Local Bubble contri- 
butions, as well as instrumental and other sources of noise. According to David Henley 
(personal communication), the halo's intrinsic intensity in the R67 band is ~ 30 counts s -1 



arcmin 2 . 



Discussion 



3.1. The Heavy Element Abundances 



Many of the calculations performed in this paper are dependent, to varying degrees, on 
the assumed abundance ratios of heavy elements to hydrogen, since we reported on observa- 
tions of atomic transitions of either oxygen (in the case of O VI emission or absorption) or an 
ensemble that includes many other heavy elements (the main source of soft x-ray emission). 
Beyond an application to our observations, heavy element abundances also influence some 
theoretical aspects of our subject matter, such as all processes that depend the rate of radia- 
tive cooli ng of the gas. While we recognize that the abundances can vary with galactocentric 



distance ( Peimbert 



1999; Martin-Hernandez et al 2002: Daflon fc Cunha 2004: Esteban et al 



1993 



2005) , and even from one place to the next at a given radius from the center (jEdvardsson et al 



Rolleston. Dufton fc Fitzsimmonslll993l ). we adopt the simplest interpretation that the 



abundances agree with the solar values. Even here, however, there are choices to be made. 
While the outcomes for determinations of the solar a bundances of elements he avier than oxy- 
gen have been reasonably stable through the years (IGrevesse fc Sauvallll998l ). very recently 
there have been some substantial downward revisions for C, N and O, based on interpre- 
tations of line st rengths in the context of detailed models of line formation in a convec- 



tive a tmosphere (IHolwegerl l200ll ; lAllende Prieto. Lambert, fc Asplundl 120021 : lAsplund et al. 



20041 ) . Of particular relevance to our work is the change in the solar abundance of oxygen 
relative to hydrogen, which has recently declined by 0.27 dex. While there has been some 
independent support for the new oxygen abundance, based on abundances found for O- and 
B-type stars (jPaflon et alll2003l ). it is discordant with models of the sound speed inside the 
Sun and the depth of its convective zone, based on the interpretations of helioseismological 
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data (IBahcall fc Serenellill2004l ; iBahcall et aj|2004j ; lAntia fc Basull2005T) coupled w ith models 
that incorporated more accurate calculations of atomic opacities ( IBadnell et al.ll2005l ). The 
problem might be solved with a higher abundance of Ne to compensate for the decr eased 
abundances of C, N and O (Untia fc Basull2005l : brake fc TestalboodJcunha et all[2006h . but 
this proposal has been met with some skepticism ( Asplund et al.l 20051 ). In view of the fact 
the solar abundance of oxygen might still be questioned, in various sections of this paper 
we will discuss the consequences of adopting either the old or new values. In cases where 
we make our own calculations, we favor the new abundances, but our use of some older, 
relatively complex models discussed in §3.41 and §3.51 were based on the old abundances. 



3.2. The Physical Properties of O Vl-bearing Gas (Isothermal Case) 



To introduce the basic concepts of our analysis of some relevant physical parameters, 
we start with the simplest case where the hot gas is isothermal and at the temperature 
that corresponds to that where O VI has its maximum ion fraction when the gas is in CIE 
(3.2 x 10 5 K). The density of electrons (n e ) in the plasma that bears the O VI ions can be 
calculated from the halo's O VI column density derived in §2.5[ intrinsic doublet intensity 
derived in §2.41 and temperature (T, assumed to be the CIE temperature, 3.2 x 10 5 K, 
but the calcul ation is relative l y inse nsitive to temperature if 10 5 K > T > 10 6 K). We use 
Equation 5 in IShull fc Slavinl (119941 ): n e = (47r/ovi)/(< ov >e Afovi), where < ov > e is the 
electron-impact excitation rate coefficient. If the observed emission originates beyond an 
extincting column of N m = 1-0 x 10 20 cm -2 , then n e = 12. 5 + j?' 2 x 10~ 3 cm -3 , but if it lies 



beyond our extreme estimates, Nm = 0.5 or 2.5 x 10 
or n e = 16.3 



jo cm 2 



respectively, then n e 



j+6.8 
-.9 



x 10" 3 cm" 3 (see Table El ) 



Table [3] also lists the thermal pressure and the depth of the emission region. The 
thermal pressure, pth, is calculated from the ideal gas law. Taking the cosmic abundance 
of the elements into account results in pth/k = 1.92n e T. The depth of the emitting region, 
Al, is equal to the column density of O VI divided by the product of the electron density 
and three calculated ratios: (H/e), which equals 0.833 in a fully-ionized plasma with a 
cosmi c composition, (O/H) , for which we adopt a value 4.57 x 10~ 4 given by lAsplund et al. 
( 120041 ). and the peak value of the fractional ion concentration arising from the balance 
of collisional ionizations and various recombination processes, fo vi(^max), calculated by 
Nahar & Pradhanl (l2003h . Our result for Al = N OY1 /{n e x (H/e) x (Q /H) m x/ n vAf^)} is 
nearly a factor of 2 larger than that presente d in previous papers ( i.e. IShelton et all (120011 )). 
because the oxygen a bundance estimates of lAsplund et al. ( 2004 ) are nearly a factor of 2 
smaller than those of Grevesse k, Anders ( 19891 ). Lastly, we calculate the time required for 
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the gas to cool, if it were to cool solely by the emission of O VI resonance line photons. (See 
Table O) The timescale estimates are inversely proportional to the oxygen abundance. The 
quoted estimates ar e appr oximately twice as large as they would have been if we had used the 



Grevesse fc Anders! ( 119891 ) oxygen to hydrogen ratio. Furthermore, the quoted timescales are 



upper limits on the true cooling timescales for the present O Vl-rich gas and do not include 
cooling of nearby, hotter gas that may, in the future, evolve through an O Vl-rich stage. 
In the following subsections, we will compare these simple predictions with those of more 
comprehensive models. 



3.3. The Extended Baseline Spectrum 



In this subsection, we create a long baseline spectrum from the halo's O VI and soft 
X-ray emission. For our comparison between the 1/4 keV and O VI intensities in § I3.4[ 
we must convert the ROSAT R12 countrate to units of intensity. In order to estimate the 
conversion rate factor, we take th e X-r ay spectrum as tha t of a T = 10 6 K, C I E plas ma, as 
determined bv iBloch et all Jl986h and IPietz etaD Jl998h : iKuntz fc Snowdenl fepQoh found 
this temperature to be the dominant temperature component in their analysis of the halo's 
soft X-ray emission (the other component is T = 3 x 10 6 K). As a result, the halo's intensity 
in ~ 110 to 284 eV photons is 2.2 ± 0.4 x 10 -8 ergs s -1 cm -2 sr _1 if the emission originates 
above Nm = 0.5 x 10 20 cm -2 and 13.4 ± 2.6 x 10~ 8 ergs s -1 cm -2 sr _1 if it originates above 
A^hi = 2.0 x 10 20 cm" 2 . The intensity of the entire 1/4 keV band amounts to only 21% of the 
intensity in the O VI doublet if the halo emission originates above extincting material having 
A^hi — 0.5 x 10 20 cm -2 and 87% of the O VI doublet intensity if the halo emission originates 
above A/hi — 2.0 x 10 20 cm" 2 . As an aside, it is of interest to compare these numbers with 
the local (Local Bubble + heliospheric) X-ray to O VI ratio, which is at least 1.4 for the 
on- filament line of sight. Thus, the local region preferentially sheds energy via the soft X-ray 
emission lines while the halo preferentially sheds energy via the O VI emission lines. We 
note that interpreting the relationship between the O VI and X-ray emission rests on the 
assumption that the same regions of space are being sampled. The OVI and X-ray data are 
from the same directions, but have differing fields of view and differing optical depths due to 
the dependence of optical depth on photon energy. As a result, the assumption is not strictly 
correct. The extent to which this is an issue depends in part on the small scale structure in 
the ISM, but is anticipated to affect our results less than the other simplifications involved 
in our modeling. 

Although a collisional ionizational equilibrium spectrum at T = 10 6 K has been found 
to fit the soft X-ray data, it underpredicts the O VI intensity (even when we assume that the 
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halo emission has been absorbed by Nm = 2.0 x 10 20 cm -2 .) As the model's temperature is 
lowered, its soft X-ray emission decreases more rapidly than its O VI emission. Therefore, it 
is possible to find a single temperature spectrum that matches both the halo's O VI and the 
1/4 keV intensities. However, the assumed temperature is significantly less than 10 6 K, and 
so does not produce the soft X-ray band ratios observed across much of the high latitude 
sky. Therefore, we will move on to more complex models which are inspired by physically 
conceivable events. 



3.4. The Expected Outcome for a Halo Supernova Remnant 



We now examine the expected consequences for a line of sight in the halo that is domi- 
nated by the effects of a single SN event. Our rationale is that in most explanations for the 
hot halo gas, the gas had been heated suddenly by an energetic event in the past. Irrespective 
of whether the energetic events occurred above the Galactic disk (i.e. extraplanar supernova 
explosions and collisions between infalling clouds and the Milky Way) or in the disk and the 
hot gas later rose into the halo, we expect the plasma to evolve from a recently shock heated 
and underionized state to a tepid and overionized state; and we expect the observationally 
determined O VI to soft X-ray ratio to be a useful diagnostic. In the following halo SNR 
simulations, we found that the ratio of O VI to 1/4 keV emission rises almost monotonically 
throughout the remnant's life, making it a diagnostic of the remnant's maturity. 

We calculated the ratio of O VI to ROSAT 1/4 keV intensities from a simulated, ex- 
traplanar supernova remnant. For the simulation, we set the ambient density, explosion 
energy, and ambient nont hermal pressur e to 0.01 atoms cm -3 , 0.5 x 10 51 ergs, and 1800 K 
cm -3 , respectively. As in ISheltonl (119991 ). the simulation included thermal conduction and 
non-CIE radiative cooling. The io nization and recombination rates were calculated using the 
Gaetz. Edgar. &: Chevalier! (119881 ) tables and the spec tra were calculated f rom t he plasma's 
non equilibri um ionization level popul ations and the iRaymond fc Smith! (119771 ) algorithm. 
We used the iGrevesse fc Anders! (119891 ) abundance tables, in w hich the oxygen to hyd rogen 
ratio is approximately a factor of 2 greater then that found by lAsplund et al.1 (120041 ) . The 
predicted O VI intensity roughly scales with the adopted oxygen abundance, while the 1/4 
keV soft X-ray spectrum is mostly unaffected. 

Here, we describe the plasma's evolution and the consequent O VI to X-ray ratio's 
evolution. When the halo supernova remnant is very young (age < 10,000 yrs), its most 
emissive portion is the hot, dense gas immediately behind the shock front. Due to the 
rapid ionizations, the recently shock-heated gas in this zone contains O VI and higher ions. 
Although this plasma emits both O VI resonance line and soft X-ray photons, its temperature 
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(T > 10 7 K) is far too high for optimum O VI resonance line emission. As a result, the ratio 
of O VI to soft X-ray intensities is less than the observed ratio (see Figure [3]). As time 
passes, the shock front slows, thus heating the gas in and just behind the shock to lesser 
temperatures and leading to a slight increase in the O VI to X-ray ratio. 

Before the remnant is ~ 50, 000 years old, both the hot dense gas near the shock 
front and the hot rarefied gas in its interior are drastically out of CIE. When the remnant 
is between 50,000 and 100,000 years old, most of its plasma comes into CIE, but not in 
the usual manner. Rather than maturing sufficiently for its ionization levels to come into 
equilibrium with the gas temper ature, the gas t emperature drops sufficiently as to match the 



ionization levels (see Figure 4 of ISheltonl (119991 )). At this time, the most important gas, that 



immediately behind the shock, is close to, but not yet in CIE. During or just before this era, 
the simulation's O VI to 1/4 keV ratio crosses that of the observed ratio (1.1, assuming that 
the intrinsic emission had been absorbed by a 2.0 x 10 20 cm -2 column, to 4.7, assuming that 
the intrinsic emission had been absorbed by a 0.5 x 10 20 cm -2 column). These observationally 
determined ratios are marked on Figure [3] by horizontal dashed and dot-dashed lines which 
cross the SNR ratio curve when the remnant is ~40,000 and ~70,000 years old, respectively. 

As the remnant continues to cool by adiabatic expansion and radiative cooling, its atoms 
begin to recombine. When the remnant is between 100,000 and 250,000 years, its shock front 
becomes too cool to produce much O VI or soft X-ray emission. Thus the shock front no 
longer outshines the bubble's interior. The remnant's O VI to ROSAT 1/4 keV ratio, now 
primarily due to the remnant's "overionized" interior, rises. Henceforth, recombinations 
from O VIII to O VII to O VI provide the remnant's center with O VI ions while reducing its 
supply of higher, more X-ray emissive, ions. In time, the temperature drops from T ~ 10 6 K 
down to T ~ 10 5 K, causing the soft X-ray emission function to plummet (although the 
O VI resonance line emission function remains nearly constant.) As a result, the O VI versus 
1/4 keV ratio continues to rise. In its final million years, the remnant is tepid (T = several 
xlO 4 K), contains some O VI ions, but few higher ions. Thus it produces an enormous O VI 
to 1/4 keV ratio, as shown in the plot. 

We observe similar trends in simulated supernova remnants having greater ambient den- 
sities, ambient nonthermal pressures, and/or explosion energies, though with some variation 
in the age when the model matches the observational ratio. Assuming that the observed 
gas can be compared to that in an undisturbed, extraplanar supernova remnant bubble, the 
time since heating and the lifetime of this gas are on the order of 10 4 — 10 5 and 10 7 years, 
respectively. The cooling timescale exceeds that calculated directly from the O VI data 
(Table [3]) because the SNR contains a reservoir of hotter, more highly ionized gas that will 
eventually transition through the O VI level. Assuming that the other possible sources of 



hot gas behave similarly to simulated halo SNRs, we suggest that the resulting structure is 
middle aged. 



3.5. The Volume Distribution Function of Temperatures 



3.5.1. The Basic Assumptions 



We move on to a more generalized picture and propose that the hot gas in the halo of 
our Galaxy is a heterogeneous mixture of regions having plasmas at different temperatures, 
created possibly by many SN events whose influences on the halo medium have merged 
together. A convenient characterization of the temperature distribution function over volume 
can take on the form of a power-law, one that extends from 10 5 K (above which appreciable 
amounts of O VI are expected in CIE) to a sharp cutoff at some high temperature T cut . Within 
the paradigm that the hot gas is created by supernova shock waves that are produced in the 
halo or that escape from the disk, the origin of T cut can be interpreted to arise from either 
a generalized limit on the supernova shock velocities or, alternatively, from our recognition 
that gases with temperatures that are too high may escape very rapidly in the form of a 
very low density, high velocity wind, which is difficult to detect. 

Since all of the observable effects of the hot gas represent line integrals of various physical 
quantities through the Galactic halo, there is some benefit in our starting with a formulation 
that describes how the many differential length segments dl within the population of discrete, 
homogeneous gas regions are distributed over temperature. We do this by specifying a 
transformation between dl and cilnT, 



where B is a distance scale factor. In a broad, statistical sense, this distribution function 
describes how temperatures are weighted according to volume fractions, but it says nothing 
about the internal electron densities n e within the length segments. 

In the development of our interpretation, we impose a simplifying constraint that the 
Galactic halo is approximately isobaric. We use this universal pressure constraint to tie 
the electron density to temperature and thermal pressure, i.e., n e = p th /(1.92kT). While 
we must accept the reality that some pressure variations can exist from one location to 
another, we can assume that the magnitude of such variations are small compared to the 
vastly different temperatures that we adopt in our model. Our representative pressure p t h 
is a free parameter that we will determine from the ratio of Iovi/Nqyi after we have solved 



dl 



BT^dlnT for T < T c 
for T > T cut , 



cut 



(1) 
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for the power-law coefficient (3 and temperature limit T cut . 



3.5.2. How the Observations Relate to the Model 



The column density of O VI is given by the expression 



iVovi = {^j (^j J fovi(T)n e dl (2) 

where (O/H) , (H/e), and fovi(T) were first used in §3.21 We can rewrite Eq. [2]in terms of 
an integral over InT of the temperature distribution we adopted, 

by making a substitution of the terms in Eq. [I] for dl, using the relation n e = p t h/ (1.92/cT), 
and noting that O VI is rare at temperatures below 10 5 K. 

We can develop similar equations for the expected strengths of the line radiation from 
O VI or the intensities of soft X-rays. For the emission from both members of the O VI 
doublet, we anticipate that 

iovi = J r vi(T)n 2 e dl , (4) 
where the emission rate coefficient roYi(T) per unit emission measure is given by 

rovi(T) = 3.09 x 10 18 (^j ^ /ovi < ^ ~ P hotons cm-V^r-^cm" 6 pc)- 1 . (5) 

As we did in §3.21 we adopt an analytical formulation of < av > (T) / (4?r) that was spec- 
ified by Shull & Slavin (1994). Again after substituting the expression in Eq. [T]for dl and 
Pth/(l-92kT) for n e we find that 

(1.92A;) 2 y lnl0 5K 

A similar equation can be expressed for the soft X-ray emission, 

2 d /-In T CU ( 

Rl2 = (£kv / r 12 [N m ,T]T^ 2 dlnT , (7) 

where ri2[Nm,T] is the emission rate coefficient as a function of temperature for 0.25 keV 
X-rays, matched to the responses of the ROSAT bands 1 and 2 and expressed in the units 
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counts cm 2 s 1 arcmin 2 (cm 6 pc)" - 1 (Snowden et al. 1997) Jd This coefficient includes the 
effect of an energy-dependent foreground absorption represented by Nm. 

3. 5. 3. Evaluation of Parameters 

To evaluate the free parameters (3, T cut and Pth/k, we compare three ratios of observable 
quantities with their expectations within our formalism. First, we may constrain the value 
of f3 by matching the observed ratio of Jovi to R\ 2 to the expression 

iovi = f::^r ovl (T)T^d\nT 
R ^ S^r 12 [N m ,T]T?-H\nT ' 

While the result depends on an adopted value of T cut , for reasonably high values of this 
quantity the effect is small. We find that for Jovi = 4680 photons cm _2 s _1 sr _1 and R\ 2 = 
0.000733 counts cm^ 2 s _1 arcmin -2 we obtain Iovi/R\ 2 — 6.4 x 10 6 (in the same units) with 
an uncertainty of 24% using the errors stated in §2J The first row of Table H] lists the derived 
values of (3 for the best and limiting values of foreground H I absorption. The top row of 
panels in Fig. @] shows how Iq, V1/-R12 changes with (3 and demonstrates that the dependence 
on T cut is not important for (3 < 1.5 and is a small effect for (3 < 2. 

A strong response to T cut arises from the ratio of high energy X-ray emission to the 
intensity at lower energies, as indicated in the second row of panels in Fig. HI where we have 
made use of an equation that is identical in form to Eq. [8] except that it compares the ratio 
of X-ray emission at 1.5 keV (the sum of intensities recorded in ROSAT Bands 6 and 7) to 
that at 0.25 keV (the sum in ROSAT Bands 1 and 2), 

R« = j;:^r 67 [N Hh T]T^dlnT 
R ^ f^r 12 [Nv h T}TP-mnT 

From the observed value of R 67 /Ri 2 = 0.04 ( §2.6[) . we arrive at values of logT cut given in the 
second row of Table HI 



lr The cosmic abundances adopted by Snowden et al. (1997) did not incorporate the recent downward 
revisions of the solar photospheric abundances of C (Allende Prieto et al. 2002), N (Holweger 2001) or O 
(Asplund et al. 2004). The correction to allow for these abundance changes should be rather small, since 
the line emission over the energy range of interest is dominated by lines from other, much heavier elements, 
such as Si, S, Mg and Fe; see, e. g., Kato (1976). 



-21 - 



< v 2 > 



Finally, we solve for pth/k by matching the observed ratio of Jovi to iVovi to the formula 

^OVI Pth 

^ovi 1.92* /^-/ovi^-^lnT 

Values of p^, found for Iovi/^ovi = 2.24 x 1CT 11 photons s -1 sr -1 with an uncertainty of 
about 50%, are listed in the third row of Table HI From the flatness of the curves shown in 
the third row of panels in Fig. HI the value of p t h is not strongly influenced by the choice of 
(3 (and has virtually no dependence on T cut ). The dominant uncertainty in p t h arises from 
the expected error in Iovi/Nqyi- 

It is important to check that the computed width for thermal Doppler broadening of 
O VI in the model does not exceed the observed widths of the actual absorption profiles 
observed in the Galactic halo. Our expectation is that 

_ / n(OVI) < v(OVl) 2 >dl _ k j^g / OVI (T)T^lnT 

fn(OVl)dl 16m p J^f ovl (T)T^dlnT ' 

where m p is the mass of a nucleon. For our derived values of (3, we find that the Doppler 

i 

parameter b = (2 < v 2 >)2 ranges from 19 to 23 km s . These values are smaller than 
the representative observed values, which are of order 60 km s -1 . Evidently, even with our 
allowance for some O VI arising from temperatures above 10 6 K, the expected profile widths 
are still much less than the kinematic broadening that takes place in the halo. As with our 
findings on pth/k, the bottom row of panels in Fig. H] shows that the outcome from Eq. [TTlis 
nearly independent of T cut . 

So far, all of our evaluations of various ratios of certain quantities have not had to 
incorporate the constant B that appears in Eq. [TJ since this term canceled out in each case. 
This constant is important, however, if we wish to relate the results from the observations 
to either a total length scale along our line of sight, L — J dl, or the values for the electron 
density, n e (T). By solving either Eqs. [31 [6] or [7] one can derive this constant, which varies 
markedly with the adopted value of (3 (as it should). Values of B for the various combinations 
of /3 and iV(H I) are listed in Table HI From an appropriate value of B we can evaluate the 
total length of hot gas along our path 

/•In Tcut 

L = B T p d\nT = B(T^ ut -10^)/(3 . (12) 

Jin 10 5 K 

Outcomes for log L (in pc) are listed in Table HI If the gas above and below the plane 
is stratified in layers parallel to the plane, then a representative vertical sight line (i.e. 
one perpendicular to the galactic disk) intersects less material than our sight line at an 
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intermediate galactic latitude of b = —47.1°. Thus, calculations for the vertical sight line 
would require that the total length scale, L in Eq. [12] be foreshortened by a factor of sin ( 1 6 1 ) . 
Similarly, iVovi in Eq. [31 Iovi m Eq. El R%2 in Eq- El together with N e and dEM(T) in the 
upcoming Eqs. [UJII31 would require reduction by the same factor. 

On the one hand, it is important to note that L represents a minimum length scale, 
since parcels of gas could be scattered over a longer distance with some unseen gas phase 
filling in the intervening gaps. On the other hand, the fraction of L that has an appreciable 
concentration of O VI is small. Figure [5] illustrates what the concentration of O VI would look 
like if all of the gas parcels at different temperatures were sequenced along a line in order 
of increasing temperature. Most of the O Vl-bearing gas is confined within a region that is 
less than about 100 pc thick, with the remaining, much greater length filled with a plasma 
that is too hot to contain much O VI. This thickness is considerably less than the 2.3 kpc 
scale height found for O VI in the Galactic halo by Savage et al. (2003). Evidently, our 
large observed ratio of Iovi/^ovi shows that the gas with temperatures of around 300,000 K 
is highly clumped, even when one acknowledges that there can be a broad distribution of 
temperatures. The internal densities n(OVI) within each clump, which in our model can 
reach as high as ~ 10~ 6 cm -3 , is about two orders of magnitude higher than an overall 
average < n(OVI) > = 1.7 x 10~ 8 cm -3 that was found in a FUSE survey of the Galactic disk 



f lBowen et all 120061 ). 



The column density of electrons along our line of sight is given by 

N e = B n e (T)T p d\nT = — P — ^^ cut )_ (13) 

y inl05K e{ ' 1.92fc (3-1 { } 

and we list outcomes for the appropriate values of (3 and Nm in Table [U 

In §3.11 we discussed the p ossibility that the abu ndance of oxygen might be higher than 



the one that we adopted f rom lAsplund et al. 



abundance value given by Grevesse fc Sauval ( 19981 ) (larger by 0.27 dex), the values of (3 



200 



4). If we were to adopt the old, higher 



listed in Table @] would increase by about 0.5 because the relative proportion of the plasma 
at temperatures that emit O VI line radiation is reduced, when compared to the higher 
temperature material that emits soft X-rays. (Recall that the predicted X-ray emission is 
not strongly influenced by changes in the abundance of oxygen, as we stated in an earlier 
footnote.) For the three values of foreground iV(H I), (0.5, 1.0, 2.0 x 10 20 cm -2 ), the old oxy- 
gen abundance changes \ogp t h/k by (+0.03, +0.03, +0.05) dex, respectively. Values of log L 
change by (+0.23, +0.09, +0.17) dex, while logiV e changes by (+0.05, +0.01, +0.03) dex. 
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3.5.4- The Total Radiative Loss Rate for the Halo 



Here, we assume that the observed line of sight provides a fair representation of the 
halo as a whole. This assumption allows us to calculate the total radiative energy loss rate 
associated with hot gas in the halo, a rate that can be usefully compared with the energy 
injection rate. 

To compute the radiative energy loss rate from the hot gas in the halo, we must know 
how the differential emission measure dEM = n 2 (T)dl varies with temperature. With our 
reformulation of dl in terms of cHnT we find that 



dEM(T) = (y^j:) 2 BT?- 2 d\nT 



(14) 



along our line of sight. The second to the last group of numbers in Table H] shows the 
logarithms of the coefficient in front of the T^~ 2 term in the equation. 

With our expression for the emission measure as a function of T, we are now prepared to 
estimate the power radiated by the gas at temperatures above 10 5 K. For a cooling function 
A N that is normaliz ed to the local product of electron and ion densities n e ni, we adopt 
the values listed by ISutherland fc Dopital (119931 ) that apply to a plasma that has a solar 
composition and nonequilibrium ion fractions in a regime of radiative, isochoric cooling, 
starting with an initial temperature of 10 7,5 K. If we were to suppose that our line of sight 
shows a fair representation of the cooling per unit area dll/dt by hot gas on both sides of 
the Galactic plane in our region of the Galaxy, we find that 



lnT cut 

d£//tft = 2sin(|6|)(0.917) / EM(T)A N d\nT 

'lnl0 5 K 



(15) 



where b is the Galactic latitude of our line of sight (— 47?1), the factor 2 accounts for 
both sides of the plane, and the factor 0.917 allows for the fact that our expression for 
EM(T) is cast in te r ms o f n 2 e whereas the normalization of the function expressed by 
Sutherland fc Dopital (119931 ) is normalized to n e rii. Values of log dU/dT are listed in the 
last group of numbers in Table HI Several important qualifications must be expressed about 
these numbers. First, we have truncated the calculation at a lower limit T = 10 5 K be- 
cause we have no ability to sense gas at lower temperatures. If an extrapolation of our 
power-law representation seems plausible, we should expect to find some additional energy 
radiated below 10 5 K. This is demonstrated in Fig. El where we have plott e d the shape of 
the integrand in Eq. [15] as a function of log T. S econd , Sutherland fc Dopital ( 1993 ) used the 
solar abundances given by lAnders &: Grevessd (119891 ) instead of the more modern values of 
Allende Prieto. Lambert. &: Asplundl (120021 ) that we adopted. Thus, if the function A N were 
to be recalculated using the newer abundances, we would find a somewhat lower cooling rate, 
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but by a factor that is less severe than the changes in the C, N and O abundances, since other 
elements are also important coolants. If indeed the older abundances are correct, the ripple 
effect from the changes in the parameters derived in §3.51 above could reduce the emission 
power by up to -0.4 dex. Aside from these possible offsets, our determination should be free 
of bias, and indeed the values of O VI and X-ray emission found here are fairly typical of 
those found along other high latitude sight lines. Nevertheless, it is important to note that 
the distribution of hot gas in the Galactic halo is extremely uneven; hence our single line of 
sight does not give a very accurate representation of the average EM on either side of the 
Galactic plane. 

The rate of cooling, dU/dt, is similar to the rate of energy input from supernova ex- 
plosions and pre-supernova winds. The average massive SN progenitor star releases 1.4 to 



2.0 x 10 ergs in wind energy bef o re it e xplodes, according to calculations in iFerrierd (119981 ) 



and iLeitherer. Robert fc Drissenl (119921 ). At the Sun's galactocentric radius, 18.6 massive 
stars and 2.6 white dwarfs explode per Myr per kpc 2 cross-sectional disk area, though many 



of these stars explode above the disk (lFerrierdll998l ). If each explosion releases 10 51 ergs of 
energy, then SN and pre-SN winds inject 7.66 to 8.06 x 10 38 erg kpc -2 s _1 into the ISM. A 
comparison with the calculated energy loss rate for our nominal case (5.37lg 69 x 10 38 er § 
kpc -2 s -1 ) shows that the majority of the injected SN and pre-SN energy is later radiated 
away by hot halo gas. Since the halo cooling rate accounts for ~ 70% of the SN and pre-SN 
wind energy injection rate, little remains to power other activities, such as large scale galactic 
winds. Half of the photons emitted by the halo will travel toward the galactic midplane and 
likely be absorbed, while the other half will travel upwards. Their absorption rate depends 
on the column density of gas residing above the emitting gas, which is expected to be less 
than ~ 1 x 10 20 cm~ 20 . Thus, most of the 1032 A and shorter photons will leave the system. 



4. Summary 

We analyzed FUSE LWRS data for i = 278.7°, b = -47.1°, finding an O VI doublet 
intensity of 4710 ± 570 photons cm -2 s _1 sr _1 . Our pointing direction is only mildly ex- 
tincted, so the observed intensity included contributions from the Local Bubble and the 
Galactic halo. Only 2° from our pointing direction, the sky is heavily extincted by a fil- 
ament residing 230 pc from the Earth. A previous observation of that direction indicated 
the Local Bubble's contribution. The difference between the intensities observed on these 
two sight lines, 4680^ggg photons cm~ 2 s _1 sr _1 , can be attributed to the halo. Given the 
extinction along our line of sight due to an A%i of 0.5 to 2.0 x 10 20 cm -2 , the halo's intrinsic 
intensity is 53501^0 to 7960 -ii2o photons s" 1 cm" 2 sr" 1 . 
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We estimated the halo's O VI column density from absorption lines seen in the UV 
spectra of extragalactic objects. We averaged values toward the 4 nearest sight lines and 
then subtracted an estimate for the expected contribution from the Local Bubble. When we 
used it and our intrinsic O VI intensity range and treated the O Vl-bearing gas as if it were 
isothermal, we found that the electron density and thermal pressure in the halo's O Vl-rich 
gas are 0.01 to 0.02 cm~ 3 and 7000 to 10,000 K cm~ 3 , respectively. 

By performing a similar shadowing analysis with the ROSAT 1/4 keV data, we deter- 
mined the 1/4 keV count rate attributable to the Galactic halo along I = 278.7°, b = —47.1°. 
It is 732 ± 142 xl0~ 6 counts s" 1 arcmin -2 , before accounting for line of sight extinction. 
After accounting for extinction, the intrinsic R12 countrate is 779 ±151 to 4770 ±930 counts 
s -1 arcmin -2 . The O VI vs 1/4 keV intensity ratio is 4.7 to 1.1, depending on whether the 
emitting gas is beyond iV HI = 0.5 or 2.0 x 10 20 cm -2 , respectively. Thus, more energy leaves 
the system through the O VI resonance lines than through the ROSAT 1/4 keV bandpass. 
In contrast, the opposite is true of the local region (Local Bubble ± heliospheric) , where 
roughly twice as much energy (at least) is radiated by the 1/4 keV emission lines than by 
the O VI resonance line doublet. 

In order to estimate the maturity of the emitting plasma, we compared the O VI to soft 
X-ray ratio with predictions for a simulated supernova remnant at various times in its life. 
Around the time that the thermal temperature throughout most of the remnant approached 
the plasma's collisional ionizational equilibrium temperature, the spectrum produced a 4.7 
to 1 ratio, coincident with the observational ratio assuming minimal extinction. Earlier in its 
life, the simulated remnant had produced a 1.1 to 1 ratio, coincident with the observational 
ratio assuming mild extinction. Specifically, these ratios were achieved when the SNR was 
70,000 and 40,000 years old, respectively. We suggest that other possible hot gas structures 
would be similarly mature when they, too, produce a similar O VI to soft X-ray intensity 
ratio. 

Using the O VI to soft X-ray ratio, we were also able to parameterize the hot halo plasma 
with a power law temperature distribution (dl = BT^d InT for 10 5 K < T < T cut ). Given the 
nominal estimates of foreground absorption (N m = 1.0 x 10 20 cm -2 ) and O VI column density 
(A"ovi = 2.09 x 10 14 cm" 2 ) along our line of sight, we found f3 = 1.48 ± 0.18, T cut = 10 6 6 K, 
and B = io _616±0 - 37 K _/3 pc. Hot gas with a temperature between 10 5 K and T cut occupies 
J dl = io 3 ' 44±0 - 37 pc along our line of sight, but the O VI- rich gas occupies a small fraction of 
this length. Assuming that our line of sight is typical of high latitude sight lines, we found 
the cooling rate for the halo (both sides of the plane beyond the Local Bubble) per unit cross 
sectional area to be dll/dt = io 38 - 73±0 06 erg s _1 kpc~ 2 . At the Sun's galactocentric radius, 
the hot halo's radiative cooling accounts for ~ 70% of the energy injected into the ISM from 



-26- 



SNe and pre-SN winds in the galactic disk and halo. The remaining ~ 30% of the injected 
energy must be split between all other energy loss processes. 

We appreciate K.D. Kuntz's assistance with the DIRBE corrected IRAS data, D. Hen- 
ley's comments on X-ray observations of the Local Bubble and Galactic halo, and S. Snow- 
den's provisions of the digitized versions of the ROSA T response curves shown in Fig. 7 of 
Snowden et al. (1997). This work was funded through NASA grants NNG04GD77G and 
NNG04GD78G to the University of Georgia and grant NAG5-12519 to Princeton University. 
This paper utilized observations obtained by the NASA-CNES-CSA Far Ultraviolet Spectro- 
scopic Explorer (FUSE) mission operated by Johns Hopkins University, supported by NASA 
contract NAS5-32985. 
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Tablc 1. Observed Intensities and 1 a Statistical Uncertainties 





Night Only 


Night Only 


Day+Night 


Day+Night 




Method #1 


Method #2 


Method #1 


Method #2 




(ph s" 1 cm" 2 sr" 1 ) 


(ph s" 1 cm" 2 sr" 1 ) 


(ph s -1 cm" 2 sr" 1 ) 


(ph s" 1 cm" 2 sr" 1 ) 


VI 1032 A 


3270 ± 460 


3150 ± 540 


3150 ± 350 


3680 ± 420 


VI 1038 A 


1440 ± 380 


2600 ± 630 


1490 ± 280 


1620 ± 380 


C II* 1037 A 


1550 ± 330 


1550 ± 400 


1700 ± 250 


1910 ± 300 


1031 A feature 






1770 ± 260 


1770 ± 270 



Table 2. Velocity with Respect to the Local Standard of Rest 





Night Only 
Method #1 
(km sec" 1 ) 


Night Only 
Method #2 
(km sec" 1 ) 


Day+Night 
Method #1 
(km sec" 1 ) 


Day+Night 
Method #2 
(km sec" 1 ) 


O VI 1032 A 


28 


25 ± 10 


22 


14 ± 10 


O VI 1038 A 


-20 


-24 ± 30 


-5 


-2 ± 10 


C II* 1037 A 


-24 


-31 ± 10 


-17 


-13 ± 10 



Table 3. Physical Parameters of Halo O Vl-rich Gas, Isothermal Case 



Assumed N u Transmission 



(10 20 cm" 2 ) 



'ovi 



N. 



OVI 



(ph s cm 2 sr 1 ) (ions cm 2 ) 



Al 

(cm -3 ) 

(PC) 



Pth/k 

tcool 

(K cm" 3 ) 
(years) 



0.5 
1.0 
2.0 



77% 
59% 



53501?*} 



6iiol38 



7960!?™ 



2.09 ± 0.84 x 10 14 



conoid 

69.8i*j;| 
012 c i + ° °° 52 

U.UlZO_g 0053 

+49.7 
-49.5 
+0.0068 
0.0069 



61.1 
0.0163 



46.8l 3 ?; 9 



6740l 2 l° 
7.36! 3 ^ x 10° 



7690l 322 ° 



6.45 



+2.74 
-2.70 



10040; 



x io € 

-4200 
-4260 



4.94l 2 ;^° x 10 6 



Note. - 



Table is "line-wrapped"; Al and t coo i appear below n e and pth/k. 
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Table 4. Parameters for the Temperature Power Law 



Foreground N(R I) (10 20 cm~ 2 ) 



Parameter 21 0.5 1.0 2.0 Eq. nr. 







1.15 ±0.20 


1.48 ±0.18 


1.95 ±0.17 


8 


loe;T r „t (K) 




6.9 


6.6 


6.4 


9 


1 ( \(V T) / 7 1 1c ( r'TTl \C \ 


IP+J 


q or _i_ n i q 


Q QQ _|_ f) 1 Q 


A no _i_ n i q 






IP J 


o.oj i u.iy 


Q no _|_ n 1 Q 


a nfi 4- n iq 

4:.uu in u. ±y 






f/3-1 


3.84 ± 0.19 


3.91 ± 0.19 


4.05 ± 0.19 




/'Irm^e 


( R4-\ 


90S 


99^ 


9^0 

Loo 


rm 




IP J 




91 A 


9 /1 9 
Z4fcZ 






IP - J 


1 Q1 


9n^ 


99Q 




log 5 (K -/J pc) 


(/?+) 


-5.37 ±0.37 


-7.12 ±0.37 


-9.95 ±0.37 


3 




09) 


-4.25 ±0.37 


-6.16 ±0.37 


-8.92 ±0.37 






(/?-) 


-3.14 ±0.37 


-5.15 ±0.37 


-7.96 ±0.37 




logL (pc) 


(P+) 


3.81 ±0.37 


3.56 ±0.37 


3.36 ±0.37 


H21 




09) 


3.62 ±0.37 


3.44 ± 0.37 


3.26 ±0.37 






(/3-) 


3.43 ±0.37 


3.31 ±0.37 


3.18 ±0.37 




logiV e (cm -2 ) 


(p+) 


19.46 ±0.19 


19.46 ±0.19 


19.50 ±0.19 


\m 




09) 


19.34 ±0.19 


19.37 ±0.19 


19.43 ±0.19 






(/?-) 


19.25 ±0.19 


19.30 ±0.19 


19.36 ±0.19 




\og[dEM{T)T 2 ~P /d In T] 


(p+) 


1.78 ± 0.06 


0.18 ±0.06 


-2.36 ± 0.06 


M 


(K /3 - 2 cm- 6 pc) 


09) 


2.88 ± 0.06 


1.12 ±0.06 


-1.36 ±0.06 






(/?-) 


3.98 ± 0.06 


2.11 ±0.06 


-0.42 ± 0.06 




\og(dU/dt) (erg s^kpcT 2 ) 


09+) 


38.69 ± 0.06 


38.76 ± 0.06 


38.72 ± 0.06 


us 




09) 


38.73 ± 0.06 


38.73 ± 0.06 


38.78 ± 0.06 






09-) 


38.77 ±0.06 


38.76 ±0.06 


38.79 ± 0.06 





-33- 



Note. - - Errors appended to the listed quantities are the formal errors that arise from 
uncertainties in the measured quantities only. They do not include systematic errors caused 
by inaccuracies in atomic physics parameters, element abundances, the assumption of isobaric 
conditions, or deviations from our temperature power-law representation. 

a (/3+) indicates value at the largest limit of (3, (/3) indicates the value at the preferred 
value of (3, and (j3— ) indicates the value at the smaller limit for (3. 




Fig. 1. — The shadowing filament and nearby sky are shown in this infrared map made 
from DIRBE corrected IRAS 100 /mi data. The filament runs diagonally across the upper 
right portion of the image. During the "shadowed" observation, FUSE was pointed toward 
the most opaque portion of the filament (278.6°, —45.3°, within the upper circle). During 
the "unshadowed" observation, FUSE was pointed toward a low opacity region just off of 
the filament (278.7°, —47.1°, within the lower circle). The FUSE LWRS aperture measures 
30" x 30", thus a small fraction of the size of the overlayed circles. 
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Fig. 2. — LiFlA spectra, binned by 0.065 A and plotted relative to the corrected wavelength 
for the Local Standard of Rest (LSR) reference frame. Top panel: spectrum from the 
satellite-night portion of the data. Bottom panel: spectrum from the full dataset. The O VI 
1032 and 1038 A emission lines and the C II* emission line appear in both the night only 
and the day+night spectra. The feature at 1031.2 A in the day+night spectrum may be the 
second order diffraction of an atmospheric emission feature of He I at 515.62 A. It, like the 
Earth's airglow emission lines (marked with crossed circles), appear much brighter in the 
day+night spectrum than in the night-only spectrum. 
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Fig. 3. — The range of observationally derived O VI to 1/4 keV intensity ratios (dot-dash line 
to dashed line) are compared with the predicted ratio from an evolving supernova remnant 
(thick solid line) and the remnant's average ratio (dotted line). The SNR's VI to soft X-ray 
intensity ratio reaches 1.1 (dot-dash line) when the remnant is about 40,000 years old and 
reaches 4.7 (dashed line) when the remnant is about 70,000 years old. The observationally 
determined ratio (1.1 to 4.7) is less than the remnant's lifetime averaged intensity weighted 
O VI to 1/4 keV ratio (8.9, dotted line). Also plotted on the graph is the time period when 
the remnant evolves from the adiabatic (Sedov) phase to the radiative phase. This phase is 
bracketed by vertical lines on the plot. 
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Fig. 4. — Plots that indicate the sensitivity of various physical parameters on the values 
of observed quantities. The three columns of panels represent solutions of the equations 
assuming different values for the foreground absorption, represented by Nm'- {left) 0.5 x 
10 20 cm -2 , {middle) 1.0 xlO 20 cm -2 , and {right) 2.0xl0 20 cm" 2 cm~ 2 . In all of the plots, curves 
for different assumed values of T cut are drawn in different styles: solid curves: logT cut = 8.0; 
dashed curves: logT cut = 6.9; dash-dot curves: logT cut = 6.6; dotted curves: logT cut = 6.4. 
Top row: The dependence of (3 on the value of Jovi/ -R12 m units of O VI photons arcmin 2 R12 
counts -1 sr _1 , as defined in Eq. [SJ From the intersection of the observed value of this ratio 
(horizontal dotted line) and the curves, we obtain the best values of (3 in each case (vertical 
dotted lines dropped down to the (3 scales at the bottom). Second row: The dependence 
of T cut on Rq 7 /Ri2 as a function of (3, as defined by Eq. [9j We used these curves to define 
appropriate values of T cut in each case by locating the one that most closely intersects the 
horizontal dotted line representing the observation at the preferred value of (3. Third row: 
The expected ratio (/ovi/-^ovi)/(Pt/i/l-92/c) as a function of (3, as expressed in Eq. [IOJ The 
intersections of the dotted lines show the outcomes for the best (3 values with the thermal 
pressures listed in Table HJ Bottom row: The expected velocity variance as a function of f3, 
as defined in Eq. [TTJ 



-38- 




Distance (pc) 



Fig. 5. — The density of VI vs. distance along our line of sight for the hypothetical case 
where all of the regions are arranged in an end-to-end sequence from low to high temperature. 
The solid curve represents our solution to the power-law temperature (or length-scale) rela- 
tion for a foreground absorption equivalent to Nui = 1.0 x 10 20 cm -2 , pth/k = 8320cm _3 K, 
and /3 = 1.48. On this curve, points are labeled according to their respective values of logT. 
The other (dotted) curves represent other values of foreground absorption, as indicated (with 
Pth/k and (3 in each case set to the respective optimum values - see Table H]). 
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Fig. 6. — Radiative energy loss rates per unit area, dll/dt, as a function logT for the 
preferred values of (3 that apply to the three possible values of the foreground N(H I) (see 
Table HJ). The curves plotted represent the product of the integrand and constants given in 
Eq.dSJ 



